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INTRODUCTION 


At the Mini-workshop on stresses and strains held at Mobil Solar on 
January 23-24, 1985, we reported on our numerical results for the in plane 
stresses and the dislocation fields in silicon ribbon. We stressed that we 
had obtained a number of convergent solutions which could be relied upon to 
provide insight into the physical phenomena involved. It was emphasized at 
the meeting that the approach to the calculation for the dislocation density 
was unique and is the FIRST such calculation ever done for any problem. In 
our model the internal structure (here the dislocation density) of the 
material itself changes due to the temperature field that is imposed and to 
the related stress field. The dislocation density changes from point to 
point and the result is a prediction of the dislocation density in the final 
ribbon. This appears to us to be a valuable aid to the improvement of 
quality of the ribbon and its ability tc function as a solar cell. We have 
had no major difficulty with convergence when the dislocation density is kept 
fixed. 

In the presentation it was honestly stated that we also had attempted a 
number of cases which did not converge. We had previously assumed that this 
lack of convergence modeled some physical phenomena and that the only open 
question was precisely what phenomena was involved. In short it did not 
bother us what phenomena it signaled. Fortunately other participants did 
not see that our divergence was not purely numerical in nature. This concern 
for the meaning of the divergence has encouraged us to perform that additional 
computations described below. There is no question that the changing dis- 
location density makes it more difficult to obtain convergence from a purely 
numerical modeling point of view. When a solution converges it usually does 
so beautifully. For example in most of the ribbon, the effective stress 
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changes by 1 Pa in 10 7 Pa from one iteration to the next after 40 iterations. 
Near the melt interface the changes are somewhat larger than 1 Pa but they are 
still very small . 

RESULTS 

Upon returning to Lexington from the Mini -workshop, we ran calculations 
on two thermal profiles closely related to the (old) EFG one. These calcula- 
tions were done for a 3 cm by 5 cm ribbon just because our program was al- 
ready set at that geometry. 

The old EFG thermal profile was 


T 4 = 437e _1 - 36x costtx + 1157e"'° 66x 
- 317e"' 47 sin(^ + £) 


The two new ones to be considered are 


T 5 = 437e" 1,36 


and 


T c = 437e" 1,36x + 1 1 57e”' 066x 
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These new profiles are shown in Fig. 1 and Fig. 2. Clearly the profile 
Tg is very smooth and Tg has a "bump" that occurs in the hot part of the 
ribbon. 


The original (old) EFG profile has an additional annealing "bump" at 



larger \ than the one shown in Fig. 1. 

After 43 iterations absolutely beautifully convergent solution for both 

new profiles are obtained while the one for the (old) EFG profile (eg. T^) 

does not converge. The dislocation density at the melt interface in these 

2 

calculations was chosen as 0.001/cm . The main results are shown in Figs. 

2 

3 and 4. The final dislocation density for profile Tg is only 7.1/cm while 
2 

it is 2375/cm for the thermal field Tg. These results illustrate the role 

of "wiggles" in the thermal profiles and their adverse effect on generating 

dislocations. Subsequently we have tried additional calculations for profile 

2 

Tg and a melt interface dislocation density of 0.0015/cm . This case also 

converges beautifully and results in a final dislocation density (at y = .9 cm) 

2 2 
of 3560/cm which is 1.5 times the value for the case when No = .001 /cm . 

IJe note in passing that there is divergence for profile Tg and a melt 

2 

interface density No = 0.01/cm . 

The "effective stress" is shown in Fig. 5. It is clear that the mechanism 
which keeps the dislocations from getting out of hand in these calculations is 
that the stress falls below the "back stress" and hence the dislocation den- 
sity does not further change with increasing x. 

A detailed look at the results also shows that the dislocation density 

2 2 

changes (almost) abruptly from No = 0.001/cm at x = 0 to N = 1.5/cm at 
x = .25 cm from the mel : interface. Perhaps this value at x = .25 cm is a 
better physical parameter to relate to ones experiences. 

For the case when the thermal profile is Tg and the ribbon is 3 cm x 

g 

5 cm, the elastic stress at the melt interface is a = -.9607 x 10 Pa 

2 

at y = 0.90 cm when the rielt interface dislocation density is No = O.COl/cm , 
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this same point has a = -.3469 x 10 Pa. When the melt interface 
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dislocation density is No = 0.0015/cm , this same point has = -.8462 x 10 

2 

Pa. The dislocation density at y = 0.9 cm and x = 2.375 cm is N = 3517/cm 

2 

for the case where Mo = 0.0015/cm . In all cases the a stresses in the 
plastic and elastic cases are already much closer to being the same by 
x = 0.25 cm. However the a stresses at x = 5 cm are actually larger in the 
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plastic domain than when considered as being elastic. They are however 
smaller than at the melt interface. 

In view of the magnitude of the dislocations generated (as shown in Fig. 
3) and the corresponding extremely small plastic strains they represent, we 
do not see how much plastic deformation can be accomodated and still produce 
good solar cells. We therefore consider the case of a small dislocation den- 
sity at the melt interface is really the one that is needed. 

This investigator believes that his ability to return from the meeting 
and at once select a different profile that converges so well, reinforces the 
notion that he has a good understanding of what really causes the divergences 
to occur. 

Fig. 6 shows the structure where a single crystal was used as the seed 
in an EFG furnace in 1974. One can note a short region of new ribbon single 
crystal and its subsequent break down into a polycrystal. This photo was 
first used in a paper by Leipold, Stirn, Fouterdyk and DeAngelis in the 
Eleventh IEEE Photovoltoris Specialists conference - 1975. It appears there 
but reproduces poony hence is included herein. 

CONCLUSIONS 

In summary we believe that results such as those shown in Figs. 3 and 4 
are of importance to this program for producing elastically good solar cells. 
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